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Abstract. 

We present an overview of recent developments in species-imbalanced ( "polarized" ) 
Feshbach-resonant Fermi gases. We summarize the current status of thermodynamics 
of these systems in terms of a phase diagram as a function of the Feshbach 
resonance detuning, polarization and temperature. We review instabilities of the s- 
wave superfluidity across the BEC-BCS crossover to phase separation, FFLO states, 
polarized molecular superfluidity and the normal state, driven by the species imbalance. 
We discuss different models and approximations of this system and compare their 
predictions to current experiments. 
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1. Introduction 

1.1. Background and motivation 

Launched by the laboratory achievement of Bose-Einstein condensation (BEC) in dilute 
alkali gases [H [2], atomic physics has undergone nothing short of a revolution. It 
has transitioned from a few-body science to the study of many-body physics featuring 
thermodynamic phases and phase transitions. Some of the most exciting experimental 
developments have included observations of lattices of quantized vortices in rotating 
BEC gases [31 IH El [6], the Mott insulator to superfiuid transition of bosons in an optical 
lattice [3 m |9], the creation of a Fermi sea of cold fermionic atoms ^0\, and the 
observation of paired fermionic superfiuids that mimic superconductivity in a metal 
or superfluidity of ^i/e [m [121 IHl [H [151 lEl [HI HE] . 

After an early period exploring weakly interacting BEC's and Fermi gases, 
recently much attention has focused on strong interactions in degenerate atomic gases, 
which can be realized, for example, by confining the atoms with a periodic optical 
potential [19], simulating a solid state crystal. A complementary way to explore strong 
interactions [201 1211 ES] in atomic gases is through the use of Feshbach resonances 
(FR) [23] , abundant in alkali atoms [2H |25l [26] , which allow the experimental control 
of the inter- atomic interactions. This crucial feature of a FR, its tunability, arises 
through the Zeeman shift (detuning, 6) of a diatomic molecular "closed" -channel 
bound state relative to the "open" -channel atomic continuum. This allows a degree 
of control (even in real time) over the strength and sign of the open-channel atomic 
interactions [271 [281 |29l [30] , that is unprecedented in other (e.g., solid-state) contexts. 

Feshbach resonance tunability has led to laboratory realizations and theoretical 
proposals for a broad range of quantum many-body phenomena. The most notable 
experimental realization is of fermion-paired s-wave resonant superfluidity, that exhibits 
the BEC-to-Bardeen-Cooper-Schrieffer (BEC-BCS) crossover between a coordinate- 
space-paired condensate of dilute diatomic molecules (for large negative detuning) to 
a Fermi-surface momentum-paired BCS regime of strongly overlapping Cooper pairs 
(for large positive detuning) [3ll[32l[33l[32[35l[36l[271[371[38l[39]. This remarkable 
controllability of the atomic interactions comes along with other tunable properties, 
including the background potential (that can be used to confine an atomic gas to one [IQ] 
or two dimensions [H], or to impose a periodic optical potential), and the atom numbers 
of the trapped atomic species. 

The focus of this review is an extremely successful set of experiments studying 
species- imbalanced ("polarized" [12]) mixtures of two hyperfine states of Feshbach- 
resonantly interacting fermionic gases. The two hyperfine states form a pseudo-spin 
1/2 system, that is thus closely related to electrons in a metal [43] or high-density quark 
matter [S], but with the advantage of the tunability of the interactions and trapping 
potential. 

Thus this research effort elucidates the response of interacting fermions to the 
frustration introduced by species imbalance, that despite attractive interactions prevents 
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a complete pairing and drives phase transitions to a variety of interesting quantum 
ground states and thermodynamic phases. A sufficiently large species imbalance 
is predicted (and observed) to destroy superfluid order of resonantly interacting 
fermions even at zero temperature. Understanding the resulting strongly interacting 
nonsuperfluid ground state is of considerable interest and may shed light on other 
strongly correlated states lacking broken symmetry such as the "pseudogap" and 
marginal- Fermi liquid regimes observed in high T,, cuprate superconductors |l5l B6] . 
heavy- fermion and Kondo-lattice systems [TTl HH] , and the phenomenology near quantum 
critical points [19] . 

By identifying the species imbalance, AA^, with the magnetization M (= AN/V 
with V the system volume) and the two-species chemical potential difference, A/i = 
fi^ — fil = 2h with an effective Zeeman energy h, this research makes a connection to a 
large body of work on a related condensed-matter system, namely superconductors under 
an applied magnetic field [13]. Thus, imbalanced Feshbach-resonant atomic systems 
allow a natural extension of superconductivity in a Zeeman field (previously limited 
to the BCS regime [50| |5T] ) across the full BEC-BCS crossover. Furthermore, unlike 
charged solid-state superconductors, where Zeeman and orbital fields are necessarily 
simultaneously induced by the applied magnetic field, atomic systems are advantageous 
as they permit independent control of these two effects, with the former created by 
species chemical potential difference, and the latter induced by the rotation of the 
atomic cloud [52ll53l[5i]. 

Experiments on polarized Fermi gases have been led by the MIT and Rice 
groups [551 ESI [571 EH] , that early on observed imbalance-induced phase-separation and 
explored theoretical predictions (such as the phase diagram [59]) for the BEC-BCS 
crossover at a finite species imbalance [601 ED E2l ESI EH ESI ESI E3 • such experiments, 
the system is typically prepared with all the atoms in one state. By applying a radio- 
frequency (RF) pulse [551 ES] , a controllable number of atoms can be transferred to the 
other state, creating an imbalanced Fermi gas. The error on the resulting imbalance 
(measured at the end of the experiment) is typically a few percent, allowing an accurate 
study of the phases of strongly interacting imbalanced fermion gases. 

Inspired by the rapid experimental progress, this line of research exploded with 
theoretical activity, leading to a wealth of predictions that continue to be explored 
in the laboratory. These include an extension of the T = microscopic analysis of 
the detuning- imbalance phase diagram [591 ES] to finite temperature [SHI ES [TOl 171] . 
proposals for numerous new phases and transitions, and an analysis of the strongly 
interacting normal state appearing at high imbalance [721 Ell El]- The universal 
phenomenology of a balanced [751 ESI EZl EHl ESI IHD] unitary Fermi gas has also been 
extended to a finite imbalance [HD IHS IHH] . Along with these theoretical developments, 
new and powerful experimental probes (e.g., RF spectroscopy [72], in-situ density profile 
measurements [57], and measurements of collective modes [^) have been brought to 
bear on this rich system, and in turn have stimulated further theoretical predictions. 
An overview of the considerable progress in the current understanding and of the 
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remaining open questions in these imbalanced resonant Fermi gas systems is the focus 
of this review. Aimed to be significantly briefer and less comprehensive, it complements 
significantly more detailed recent reviews by Gurarie and Radzihovsky [20] . by Bloch, 
Dalibard, and Zwerger [19], by Giorgini, Pitaevskii, and Stringari [21] and by Ketterle 
and Zwierlein [22]. 

The present review provides an overview of the current status of the field of 
imbalanced fermion gases, and is considerably broader than our recent research 
article [66] that only provided the details of our own theoretical calculations of the 
T = mean-field phase diagram. In this review, we aim to avoid the details of analyses 
while discussing the present-day understanding of the T = mean-field theory [591 EE] 
and the results of finite temperature mean-field and large-N calculations in Sec. [2l the 
nature and stability of strong fluctuations in the FFLO state in Sec. [3], the predictions 
of the local density approximation in Sec. HJ and the status of the most up to date 
experimental observations on this system in Sec. [51 In Sec. [6l we conclude with a few 
remarks about the future directions of this exciting subject in the physics of correlated 
fermions. 



1.2. Two-body Feshbach-resonant scattering 

Before turning to the analysis of states at finite density, we briefly review the Feshbach- 
resonance mediated interaction, that (at finite density) leads to superfluidity in fermionic 
atomic gases. From a broad perspective, the low-energy atomic scattering and therefore 
(through the T-matrix) the effective inter-atomic attraction is strongly influenced by 
the presence of a bound (or quasi-bound) state in the two-atom potential. More 
microscopically, Feshbach resonant scattering in the "open" channel takes place when 
there is an energetically- near by "closed" channel, that is weakly coupled to it by the 
hyperfine interaction. Crucially, due to a Zeeman sphtting between the bound state 
(closed channel) and the scattered states (open channel), an external magnetic field can 
be used to change the bound-state energy relative to the continuum of the scattering 
atoms, and thereby tune the effective interactions between scattering atoms. 

The Feshbach-resonant scattering of the open-channel fermions is characterized by 
the s-wave scattering amplitude fo{k), that by unitarity and analyticity [84J is required 
to take the form 

with as the scattering length and tq the effective range parameter that is negative in the 
resonant regime [20]. As the closed-channel bound-state energy is tuned through zero 
by adjusting an external magnetic field B, the scattering length diverges as 

as = %g(l- q^^qJ ^ (2) 

with abg the background scattering length, Bq the resonance position (at which the 
bound-state energy vanishes) and B^j the resonance width. 
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A microscopically faithful way to capture the Feshbach resonance phenomenology 
is through the so-called two-channel model [351 ESI 122], in which the closed channel 
molecules appear explicitly in the Hamiltonian and give the scattering amplitude in 
Eq. ([1]). In practice, however, present-day experiments on superfluid Fermi gases can 
only access low energies (densities), that can be described by a simpler, more universal 
one-channel model to which the two-channel model reduces in the low-energy, fc|ro| ^ 1 
limit [20| |85] . We thus focus on this experimentally relevant limit, described by the 
one-channel model Hamiltonian 

k,cr kqp 

for the two species (cr \) of open-channel fermions created by the anticommuting 
operator c|^^, with the single-particle energy = h'^k'^/2m, mass m, and volume V. 
The attractive interactions are parameterized by A < 0. 

Through a standard T-matrix scattering calculation, that gives —fo{k) = 
the pseudo-potential parameter A can be related to the experimentally determined, 
magnetic field dependent |25l |26] scattering length 

ATia^f9 ~ A^ \7^2^' ^ ' 

k 

where the ultraviolet-divergent second term is regularized by a microscopic momentum 
cutoff scale A ~ l/ii set by the closed-channel molecular extent d. This gives (with 
h = 1 throughout) 

/ 47r 2A ^ X /rx 

«s A = — + — = — Ar, 5 
\mX -K J Ait 

= ^ ^ (R) 
47rl-A/Ae' 

where A^ can be called the renormalized coupling and Ac = is the critical coupling 
A at which the scattering length diverges. The above relation allows the definition of the 
model and therefore a reexpression of physical observables in terms of the experimentally 
defined (UV-cutoff independent) scattering length Og. 



1.3. Resonant Fermi gas at finite density: limits of validity 

In contrast to the two-body problem discussed above, a resonant many-body system at 
finite density presents a formidable challenge, that in most cases can only be treated 
approximately. Although many uncontrolled (but illuminating) approximations have 
appeared in the literature [311 EHl EHl [3H] dating back to the pioneering works studying 
the BEC-BCS crossover [311 [32l [33] , such one-channel model studies across the unitary 
point can only be trusted qualitatively. The reason is that the one-channel model at 
finite density n exhibits only a single dimensionless gas parameter na^ (or, equivalently, 
/cpOs, with the Fermi momentum kp), a measure of the resonant interaction strength, 
that diverges upon approach to the Feshbach resonance, where |as| — )■ oo. Thus, a 
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perturbative expansion in kpas is precluded sufficiently close to the resonance and is 
useful only far from the resonance, deep in the BCS and BEC limits, where /cfIc^sI ^ 1- 

Nevertheless, trustworthy treatments of such strongly interacting, scale-invariant 
systems are indeed possible, taking a cue from the studies of critical phenomena [861 
[871 [88] . where analogous challenges have been surmounted almost 40 years ago. The 
approach is to "deform" the physical model to one that exhibits a mathematical small 
parameter, such as the deviation e = duc — d of spatial dimension d from an upper-critical 
dimension, due-, above which the behavior is simple, or a l/Nf expansion about a large 
number of atom flavors, Nf, with the exactly solvable Nf oo limit formalizing and 
justifying earlier mean-field approximations. There has been considerable success in 
applying these field-theoretic methods to a unitary Fermi gas[82l [831 ESI HO]- These 
methods have also been embellished in uncontrolled "self-consistent" resummation 
schemes, that do not always lead to an improved description [91]. From the related 
renormalization group point of view, in the solvable vacuum (two particle) limit of the 
previous subsection, a Feshbach resonance with |as| — > oo corresponds to a critical 
point [821 [83], that separates two distinct "phases" (vacua at zero density) on the 
BCS (positive detuning, Os < 0) and BEC (negative detuning, Os > 0) sides of the 
resonance. At a finite density n, this transition between the corresponding BCS and BEC 
paired superfiuids is converted to a smoothed crossover, analogous to a paramagnet- 
ferromagnet crossover in a magnetic field. 

Another useful approach to a controlled theory of the BEC-BCS crossover is through 
a narrow Feshbach resonance |27l [20] . described by the two-channel model [23]. In 
the latter, in addition to the gas parameter /cfc^s; there appears another dimensionless 
parameter. 



that is the dimensionless ratio of the width of the resonance, Fq (proportional to B^, 
set by the strength of the Feshbach resonance coupling, the hybridization amplitude of 
an atom-pair with a molecule) to the Fermi energy ep. This resonance- width parameter 
7 naturally allows a distinction between wide (7 ^ 1) and narrow (7 ^ 1) resonances. 
Equivalently, these are contrasted by whether, upon growth near the resonance, the 
scattering length as{B) reaches the effective range |ro| first (the broad resonance) or 
the atom spacing kp^ ~ n~^^^ first (the narrow resonance). A quantitatively-accurate 
description of the full BEC-BCS crossover, perturbative in 7, can thus be obtained 
for narrow Feshbach resonances, making them attractive from the theoretical point of 
view. The ability to treat narrow resonant systems perturbatively physically stems from 
the fact that such an interaction, although arbitrarily strong at a particular energy, is 
confined only to a narrow energy window around the resonance energy [271 [20] . 

In practice most s-wave Feshbach resonances studied to date are broad, 7^1, 
with one notable exception discussed in Ref. [92]. However, even for broad resonances 
(characteristic of a featureless attractive potential, where due to lack of a large potential 
barrier no long-lived resonant state exists at positive energy) 7 can be thought of as 




(7) 
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a mathematical expansion parameter about the solvable narrow resonance (7 0) 
limit, analogous to the and e expansions. Furthermore, because of the density 

dependence 7 ~ in principle 7 can be reduced by working at higher atomic 

densities. 

1.4. Balanced BEC-BCS superfluid crossover 

Before turning to the main subject of polarized (imbalanced) resonant Fermi gases, we 
briefly review the salient points of a balanced gas and the associated BEC-BCS crossover. 
Since the physical interaction is everywhere attractive (although the scattering length 
changes sign), the ground state of the system is expected to be an s-wave superfluid for 
the full range of detuning. The latter will transition to a nonsuperfluid thermal state 
at a sufficiently high temperature (tunable via the FR), or to a Mott or Bose-glass 
insulator when subjected to a sufficiently strong commensurate periodic potential or 
quenched disorder [8]. 

The qualitative picture for the superfluid state across the resonance follows from 
the aforementioned solvable narrow resonance limit [271120] . in which the closed channel 
molecule is long lived and can therefore be treated as an independent particle even 
for positive detuning. For detuning larger than twice the Fermi energy, closed-channel 
molecules are too energetically costly and most of the atoms are in the form of open- 
channel fermions, forming a weakly BCS-paired Fermi sea, with an exponentially small 
molecular density, induced by the weak Feshbach resonant coupling. The BEC-BCS 
crossover initiates as the detuning is lowered below an energy scale of order of 2ep, where 
a flnite density of atoms begins to bind into Bose-condensed closed-channel molecules, 
stabilized by the Pauli principle. The resulting molecular (closed-channel) superfluid 
coexists with the strongly-coupled BCS superfluid of (open-channel) Cooper pairs, that, 
while symmetry-identical and hybridized with it by the Feshbach resonant coupling, is 
physically distinct from it. This is made particularly vivid in highly anisotropic, one 
dimensional traps, where the two distinct (molecular and Cooper-pair) superfluids can 
actually decouple due to quantum fluctuations, suppressing the Feshbach coupling at 
low energies [93] . The crossover to the BEC superfluid terminates around zero detuning, 
where the conversion of open-channel atoms (forming Cooper pairs) into closed-channel 
molecules is nearly complete. In the asymptotic negative-detuning regime a true bound 
state appears in the closed-channel, leading to a positive scattering length and a two- 
body repulsion in the open-channel. In between, as the position of the Feshbach 
resonance is tuned through zero energy, the system is taken through (what would at 
zero density be) a strong unitary scattering limit, corresponding to a divergent scattering 
length, that is nevertheless quantitatively accessible in the narrow resonance limit, where 
7 ~ l/(A;i?|ro|) plays the role of a small parameter. 

Although the above features of the crossover are no longer expected to vividly 
appear in the experimentally relevant broad resonance limit, it is useful to keep 
this qualitative picture in mind. The detailed description of the broad resonance 




Figure 1. Order parameter A (panel a) and chemical potential /i (panel b) as a 
function of — (fepOs)^^, with the BEG regime on the left and the BCS regime on the 
right. The solid black lines are the mean-field results for fi and A, following Ref. [94]. 
The red circles include the 0{1/Nf) corrections [83], evaluated at Nf = 1. The star 
symbols at unitarity are the results of the quantum Monte Carlo calculation of Ref. [77) . 



emerges when the Hamiltonian is treated in the mean-field approximation (although the 
aforementioned e- and l/iVj- expansions can also be utilized to systematically correct 
it). At T = 0, a standard mean-field analysis (that is exact in the large Nf limit |82l l83] ) 
of = H — fiN gives the grand-canonical ground state energy density 

Eg= -^ + V-'J2i^k~E,), (8) 

k 

with = efc — /i and i?fc = + A^, and in the second line Eq was reexpressed in 
terms of the scattering length, thereby also leading to a finite pairing integral (the second 
term). A minimization of Eq{^!S.^\x) over the variational superfiuid order parameter A 
ensures its expression in terms of the fermion pair A = A(c4.(x)c-|-(x)), i.e., gives the 
gap equation, and the chemical potential enforces the total fermion number equation 

The integrals in the mean-field ground-state energy Eq. ([8]), and the gap and number 
equations, can be expressed in terms of elliptic integrals |9l], and the simultaneous 
numerical solution of the gap and number equations yields the mean-field A and [i as 
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a function of the detuning parameter plotted as solid lines in Fig. [TJ The red 

circles include the leading correction, at Nf 1, following Ref. [83], while the 

star symbols at the unitary point are quantum Monte Carlo results [77]. 

It is also illuminating to obtain limiting analytical expressions to the mean-field 
results. In the BCS regime (1 < —]^), /i > and A ^ /x, one finds (defining c = 
with the three-dimensional density of states A^(e) = Cy/e) [66] : 



m 



15^ \2 Se-^fiJ 



4.a,- " 

that is non-analytic in A, reflecting the well-known fact that an arbitrary weak attraction 
yields a pairing instability. To leading order in A ^ /i, this gives the standard BCS 
expressions for the gap (0 = and number {n = — ^f-) equations 

A ^ Abcs = Se'^^ exp (^rj^) , (11) 

/i ^ eF. (12) 

Here, the Fermi energy ep is defined in terms of the particle density via n = ^ce^'^. 
In the asymptotic negative-detuning EEC regime of < —1, in which ^ < and 

^ A, one instead obtains the analytic Landau form 

^ A 9 1 A 7, • ,9/ TT A^, , , 

Eg ^ -^A^ + -c^ A \n + -— , 13 
4:7ias 2 32 /i"' 

that leads to the gap and number equations (for A ^ |/i|) 

f^^ - (14) 



16 

V oirkpas 

with /i equal to one-half the two-atom binding energy in the BEC limit. 

As an aside, we note that expressions for the ground state energy, the chemical 
potential, and the order parameter (gap) are all well-defined in the unitary — )■ oo 
limit, respectively given by universal forms 

3 

Eg = 0.5906-neF, A = 0.6864eF, /i = 0.5906eF, (16) 
5 

all expressible in terms of the only energy (or equivalently, density n = kp/STc"^) scale 
eF of the problem, with the numerical prefactors approximations (to lowest order in the 
1 /A^/-expansion [821 [83]) for what are expected to be universal numbers characteristic 
of a unitary Fermi gas. 

The above analysis can be generalized to finite temperature, giving the transition 
temperature to the normal state, at unitarity also scaling with ep (but which 
becomes inaccurate in the deep-EEC limit [32] )• A finite temperature also introduces an 
additional energy scale, at unitarity leading to the generalization of the above numerical 
prefactors to dimensionless universal scaling functions of ep/T. 
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2. Imbalanced resonant Fermi gas 

2.1. Model 

With tliis background in mind the imbalance in the number of the two species is 
incorporated by introducing two distinct chemical potentials, /io- = to impose 

separately conserved fermion numbers = {N^,N^), or equivalently total fermion 
number N = + and species imbalance AA^ = A^^ — A^^^. The corresponding 
grand-canonical Hamiltonian if^^ = H — jiN — hAN is 

k,(T kqp 

with fi^ = fi + h and fi^ = fi — h related to the total-number chemical potential 
fi and species-imbalance chemical potential h (an effective Zeeman energy). The 
imbalanced resonant Fermi gas thermodynamics as a function of N, AN,T,as, i.e., the 
extension of the BEC-BCS crossover to a finite imbalance can then be computed by a 
variety of theoretical techniques, including quantum Monte Carlo [631 EH], mean-field 
theory [66l ETJ [65l EH [68] , the large- A^/ (fermion flavor) [82l[83l[90] and e expansions [89] . 
By studying Eq. (ITTll we have already specialized to the wide Feshbach resonance limit. 
However, this problem can also be studied in a controlled narrow resonance regime by 
working with a two-channel model as done in Refs. [591 EH] . 



2.2. Analysis at zero temperature 

We now outline the resulting picture (that is in qualitative agreement with 
experiments [55l[56]) using standard mean-field theory [66] in the broad-resonance limit. 
At zero temperature this amounts to a minimization of the mean-field ground-state 
energy Eq with respect to the variational parameters (defined below) at fixed and h. 
The derivatives with respect to fj, and h 

N^_dEG AI^^_dEG , . 

V dfi' V dh' ^ ' 

then respectively give the total atom number, N = + and species imbalance, 

AA^ = — N^, imposed in cold-atom experiments. Our variational ansatz assumes 
pairing of the form 

A(x) = AQe^Q•^ (19) 

with Aq and Q variational parameters (with Q 7^ to allow for the possibility of the 
simplest Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) [961 [97] typ^ state). With this, one 
can derive the ground-state energy density Eg{Aq, Q, fi,h) following Ref. 



/d^k 
j^,[E^^e{-E^+) + E^^Q{-E^^)], (21) 
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with 6(x) the Heaviside step function and where we defined 



Eh 



Q 



2 



fe + ^)' + ^Q ' ' (22) 



5m 



1/2 



E^,± = EkT{h-^^), (23) 

which can be solved for the stationary Aq and Q by supplementing the gap equation 
with the condition 

Although the mean-field approximation is expected to be quantitatively inaccurate in 
the strongly-interacting unitary regime, it does provide a consistent qualitative picture 
and a natural starting point for more accurate theoretical methods. 

An important point to keep in mind, for a proper analysis of Eg[Aq,€1, fi, h), is 
that not every simultaneous solution of the gap and number equations corresponds to a 
physical ground state of the system; the key additional criterion is that the solution Aq 
must also be a minimum of Eq at fixed /i, h. The verification that an extremum solution 
is indeed a minimum is particularly essential when (as is the case for a polarized Fermi 
gas) there is a possibility of a first-order transition, where Eq exhibits local maxima 
and saddle points that separate its local minima. By working instead directly with 
the gap equation and imposing fixed and AA^ from the start, several authors have 
been lead to qualitatively and/or quantitatively incorrect phase diagrams [651 EH 169] , 
mistakenly associating local minima, saddle-points and/or even maxima with a ground 
state, and misidentifying first-order phase boundaries with the spinodals, as discussed 
in Refs. [671 [99]. 

Therefore, in our experience the safest path to the correct phase diagram is to 
work in the grand-canonical ensemble, use the global minimum of the free energy to 
map out the phase diagram at fixed fi, h, and only then constrain it to the imposed 
total number and imbalance densities; the in-satisfiability of this final constraint is a 
signature of phase separation, resolved by an inhomogeneous coexistence of two phases 
in proportions imposed by the two density constraints. 



2.3. Phase diagram 

2.3.1. Zero temperature With these cautionary remarks in mind, as detailed in Ref. [66] 
Eg{Aq,CI, n, h) can be analyzed analytically in a number of asymptotic regimes and 
numerically throughout. The result of its minimization directly leads to the zero- 
temperature phase diagram [66] at fixed density and chemical potential difference in 
Fig. [2] and at fixed density and imbalance in Fig. [3l The BCS regime of both figures is 
shown in Fig. |H 

The qualitative features of the phase diagram can be readily understood. At 
sufficiently small Zeeman energy h the fully-paired P = superfiuid state (that has 
Q = 0, with Aq=o = a = Abec-bcs) is stable across the full BEC-BCS crossover. As 
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Figure 2. Mean-field phase diagram of imbalanced Fermi g function of 

inverse scattering length and chemical potential difference (Zeeman energy) h, showing 
regimes of an imbalanced magnetized superfluid (SFm) bounded by h„i and /ici, and 
a strongly-interacting normal Fermi liquid. The region denoted SF is the balanced 
BEC-BCS crossover superfluid phase. The central shaded region is a regime of phase 
separation, where the gas inhomogeneously coexists between two phases. The FFLO 
region is indicated by an arrow and bounded, to the right of the blue point by /ifflo 
and hc2', however it is too small to be seen on this scale. A zoom- in of this region is 
shown in Fig. |4l The red point at (—2.37, 6.89) is a quantum tricritical point (QTP) 
separating the first- and second-order transitions from the SFm to the polarized Normal 
phase. The green point at l/{kpa) « 1.01, ~ 1.5 is a critical end-point. 

first shown by Sarma [51], in the BCS regime the paired state cannot accommodate any 
polarization (except for the Q FFLO state, which we shall neglect for the moment). 
This follows because the BCS s-wave paired state involves pairing at the Fermi surface 
and is fully gapped. However, in the presence of a population imbalance, the non- 
interacting Fermi surfaces have different volumes, determined by the Fermi momenta 
fcpt ci-nd kpi, disrupting pairing. With increasing h and at fixed fj,, then, the BCS state 
undergoes a first-order transition to an unpaired magnetized normal state. At imposed 
average density (atom number A^) and fixed h, this leads to phase separation for h > hd, 
as shown in Fig. |2j Since the polarization only becomes nonzero for h > h^i in the BCS 
regime, it is clear that, at fixed polarization P, one finds phase separation at arbitrarily 
small P, as shown in Fig. |3l Thus, an addition of excess spin-'|" atoms to the paired 
BCS state leads to their spatial phase separation from the paired gas. In contrast, as 
illustrated in Fig. [3l in the deep BEC regime the imposed polarization is accommodated 
by a transition of the P = molecular superfluid to a polarized superfluid, SFm, that is a 
homogeneous ground state of condensed molecules and a single-species (fully-polarized) 
Fermi gas. 

The low-/i stability of the BEC-BCS superfluid is limited by the lower-critical 
Zeeman field, that in the weakly-paired BCS regime {—l/k-pa^ ^ 1) is of the order 
hci ~ Abcs/v^5 exponentially small in l/Zcp^s (this is the Clogston limit |5D])- We define 
hci to be the Zeeman field above which the system phase separates. With decreasing 
detuning, it grows to hd ~ ep around the unitary crossover regime (—1 < — I/ZcfOs < !)• 
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Figure 3. Mean-field phase diagram of imbalanced Fermi gases, as a function of inverse 
scattering length and imbalance P, showing regimes of magnetized (imbalanced) 
superfluid (SFm), FFLO (shaded red, bounded by Pfflo and Pc2) and normal Fermi 
liquid. The central shaded region is a regime of phase separation, and the red point 
at (—2.37, 1) is a quantum tricritical point (QTP). The FFLO phase is only stable to 
the right of the blue point. The green point at l/{kFa) « 1.01, P = is a critical 
end-point. 



Deep in the BEC regime (l/fepCs ^ —1), the stabihty of the BEC-BCS superfluid is 
instead controlled by the molecular binding energy hm ^ Ei, l/ma"^. 

In the opposite limit of large Zeeman field, h > hc2 (at fixed A^), a normal (unpaired) 
state is clearly stable, with hc2 exponentially small (but larger than hd) in the BCS limit. 
With decreasing detuning, hc2 grows through ~ ep in the crossover regime and increases 
further with the binding energy (while satisfying hm < hd < hc2) deeper in the BEC 
regime. The high field-induced depaired state is a partially (two Fermi surfaces) or 
fully (single Fermi surface) polarized normal Fermi liquid, A = 0, P 7^ 0, respectively 
depending on the Zeeman field strength h and the interaction strength 1 / kpa. Deep in 
the BEC regime this normal state is fully polarized, characterized by a single-species 
Fermi surface. 

In the deep BCS regime, the interactions are weak and the large h depaired state 
is clearly a weakly renormalized polarized Fermi liquid, while in the deep BEC regime 
the only depaired state is fully polarized (P = 1), that is, therefore, noninteracting. In 
contrast, in the crossover regime (—1 < — l/fcpOs < 1) the interactions are resonant and 
strong, and therefore cannot be treated perturbatively. Nevertheless one can explicitly 
demonstrate using many-body l/A^'j-expansion techniques, that (earlier suggestions 
based on RF spectroscopy notwithstanding [72]) this depaired state remains a normal 
Fermi liquid, albeit a strongly renormalized one pO]. For example, at unitarity and 
P — 7- 1^ it is characterized by a suppressed quasi-particle residue Z = 0.47, an enhanced 
effective mass ^ 1.8m, and a downward chemical potential shift 6fi^ = — 1.46eF. 

In the presence of a population imbalance, another possible ground state is the 
FFLO phase, in which atoms pair at a nonzero center of mass momentum Q ~ fcpt — fcpj, 
to accommodate the spin-"|" and spin- J, Fermi-surface mismatch. Within mean-field 
theory, and assuming the simplest FFLO-type state, Eq. ( |T9l) . in the BCS regime the 
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above-mentioned first-order transition with increasing h is into the FFLO phase. With 
further increase in h the FFLO phase undergoes a continuous transition into the normal 
phase yielding the phase diagram at fixed h and fixed density shown in Fig. H^, and the 
phase diagram at fixed P and fixed density shown in Fig.HjD. We observe that the FFLO 
phase is favored only for a thin sliver of h or P, and for sufficiently weak interactions 
/cfIcsI < 2 in the BCS regime. In SecjSl we discuss FFLO-type states in more detail and 
the strong effect of fluctuations on them. 

Turning to the BEC regime, as illustrated in Fig. [21 at large negative detuning 
{—1/kpas < —1) and at intermediate Zeeman fields, < h < hd, an imposed species 
imbalance can be accommodated by a homogeneous polarized (imbalanced) superfiuid 
state, SFm (that we originally dubbed "magnetized superfiuid" [59] ) p^O] . This is not 
surprising as in this regime of strong attractive interactions pairing takes place in real 
space and a nonzero imbalance corresponds to the addition of a low density of excess 
fermions to a fluid of tightly-bound molecular bosons. As discussed in Ref. |101j . this 
is analogous to the well-studied problem of ^He-^He mixtures (and, generally, boson- 
fermion mixtures |102[ 1103] ). with fermionic '^He atoms corresponding to the excess 
spin-'t" fermions and known to exhibit a homogeneous phase (miscible in ^He) below a 
concentration of about 7% (at T = 0) |104l 1105] . 

In the simplest approximation, the SFm phase is characterized by a ground state 
that is the product of a fully polarized majority Fermi sea and a diatomic molecular 
condensate. This phase is therefore quite novel, as it is a hybrid state displaying strongly- 
paired superfluidity and gapless fermionic excitations around a single polarized Fermi 
surface, and is thereby expected to display features (e.g., linear low-temperature heat 
capacity, superflow, etc.) characteristic to both. The SFm is separated from the balanced 
P = superfiuid by a continuous transition at h = hm, the latter set by the binding 
energy. 

However, as illustrated in Figs. [2] and [3l over the large shaded portion of the 
phase diagram, at intermediate detuning and Zeeman field (and polarization), it is 
not possible to accommodate a species number imbalance with a homogeneous, single 
component global minimum of Eq- This happens whenever two competing ground- 
state minima become degenerate, and are characterized by two distinct densities and 
polarizations for the same critical values of chemical potentials. Thus, as is standard 
for first-order transitions, it is not possible for the system to achieve a homogeneous 
state with intermediate values of the imposed densities [591 EE] . For the corresponding 
range of parameters in the shaded region the system thus phase separates [61| into two 
coexisting degenerate ground states with different densities in chemical equilibrium, such 
that the total imposed number and polarization constraints are satisfied. The resulting 
phase-separated state can be explicitly accounted for by generalizing the ground-state 
ansatz to include the possibility of such an inhomogeneous mixture [66] . 




One feature of the above phase diagram that has attracted recent attention is the 
quantum tricritical [106j point (QTP) [681 [66], appearing at — 1/(/cf«.s) — —2.37 in the 
mean-field approximation and denoted as a red point Fig. [31 It separates a first-order 
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(to the right) and second-order (to the left) P-increasing transition out of the polarized 
superfiuid, SFm- Upon decreasing P from unity, the QTP separates two possible fates of 
an added spin-down fermion immersed in the polarized Fermi sea: To the left of the QTP, 
the magnetized superfiuid phase appears, as added spin-down atoms form tightly-bound 
pairs with spin-up atoms, miscible with the remaining spin-up Fermi sea. In contrast, 
to the right of the QTP, the added spin- J, atoms will form molecular pairs, that are 
immiscible with the remaining majority Fermi sea (leading to phase separation to the 
right of the QTP). The third possibility for the added spin- J, atoms is to form a spin- J, 
Fermi sea. This is realized to the right in the phase diagram, beyond —-j^ — —0.1, as 
seen in Fig. [31 where the phase separation region meets the normal Fermi liquid region 
at P ^ 1. 

A similarly interesting region of the phase diagram is the point (marked green in 
Figs. [2]and|n]) at which the SFm state emerges (a critical end point, at which a line of 
second order transitions terminates at a first-order transition). This point is located on 
the P = axis, at — l/(A;Fas) ~ —1.01 in the mean-field approximation and separates 
first- and second-order transitions, with increasing P, out of the fully paired, balanced 
BEC-BCS superfiuid. 

The above outlined mean-field phase diagram is in qualitative agreement with recent 
QMC calculations [95], with the quantitative discrepancies primarily traceable to the 
neglect, within mean field theory, of normal state interactions |107] . The latter lower the 
normal state free-energy and corresponding Pc2 below their mean-field values. Indeed, 
recent experiments |108[ 1109] have observed the Pc2 phase boundary (but not the FFLO 
phase) in qualitative agreement with Fig. |3l finding, at unitarity, Pc2 — 0.3 — 0.36 (Fig. 2 
in Ref. |109] ) significantly below its mean-field prediction of P"^^ ^ 0.93, but quite close 
to the result of the large- iV expansion t83J of pI^^'^ ~ 0.302 and the QMC result [95] 
of P,*^^^^ ~ 0.4. 

However, experiments have not yet detected the Pd boundary separating the SFm 
and phase separation regimes, although the finite-temperature SFm phase has been 
observed; mean- field theory and QMC predict similar values for Pd. To compare, 
we note that the tricritical point at P = 1 occurs at {k-pa^)'^ = 2.37 in MFT and 
(/cpcts)"^ = 2.14 within QMC |95]. Additionally, the onset of the SFm phase occurs at 
(fcpas)"^ = 1.01 in MFT and at (fepas)"^ = 0.53 within QMC [95]. 

One regime in which the problem simplifies (while still exhibiting strong-correlation 
effects) is the extreme polarized limit P — t- 1^ of a small density of spin-J, atoms in a 
Fermi sea of spin-'|" atoms, with the limiting case of one "impurity" spin-j, atom immersed 
in the majority (spin-'|") Fermi sea. This problem has been studied using variational 
wavefunction |110[ 1111] and Monte Carlo |112] techniques, finding that one spin- J, Fermi 
atom interacts strongly with the spin-'|" Fermi sea (forming a "polaron"), as reflected in 
values of the quasiparticle residue Z, effective mass m^: and chemical potential Sfi^ that 
deviate from that of free fermions. These polarons have been subsequently observed 
in experiments studying the P — )■ 1~ limit, obtaining (at unitarity) Z ^ 0.39 [73] , 
~ 1.17m [74j, and ~ — 0.64eF [73], values that are in good agreement with QMC 
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Figure 4. Panel a is a zoom-in of tlie phase diagram of Fig. [2] at fixed h and inverse 
scattering length. Panel b is a zoom-in of the phase diagram of Fig. [3] at fixed P and 
inverse scattering length. The horizontal axes are the same in the two panels. In each, 
the FFLO phase is only stable to the right of the blue point (as in Figs. [2] and [S]). 



and variational theories. Upon increasing the strength of the attractive interactions 
into the BEC regime, around 1/kpa ^ 1.27 a phase transition is predicted [113j from a 
polaron-hke state of the minority atom dressed by the majority Fermi sea (characterized 
by a nonzero fermionic pole with a finite residue Z < 1) to a bosonic molecule dressed 
by the (majority) polarized Fermi sea and characterized by a vanishing residue of the 
fermionic pole (Z — > 0). The properties of the latter state have been recently analyzed 
by a molecular variational wavefunction [113(1114] . thereby constructing a phase diagram 
analogous to Fig. [21 

In addition to the basic (and well accepted) qualitative features of the imbalanced 
resonant paired superfluid, summarized by Figs. [2|3| and [11 a number of additional 
phases have been suggested that are not captured by the simplest mean-field ansatz. 
These include breached pair phases [60] , deformed Fermi surface (nematic) phases [115[ 
I116j . possible further instabilities of the fully-polarized Fermi sea of the SFyi phase [lOlj . 
and a number of more exotic fluctuation-driven FFLO-type phases that we discuss in 
Sec. [31 



2.3.2. Non-zero temperature Soon after the zero-temperature phase diagram in Fig. [3] 
was established [591 [66], the analysis was generalized to a finite temperature [68 ] [69 ] [70 ] 
[71] . as illustrated in FigJSl As usual, the finite-temperature thermodynamics and the 
corresponding phase diagram emerges from the minimization of the free energy F, that 
for the imbalanced resonant Fermi gas (because of a number of possible extrema) must 
be done with care similar to that of the zero-temperature case. While the T = phase 
diagram is largely established (despite an ongoing search for the FFLO phase, remaining 
questions about its extent to smaller imbalance, and the need for more detailed studies 
of the SFm), the status of the finite-T phase diagram in the strongly-interacting regime 
is more uncertain. The basic reasons for this can be traced to the original analysis by 
Nozieres and Schmitt-Rink [32] on the finite temperature balanced paired Fermi gas and 
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Figure 5. A schematic phase diagram for fixed detuning, illustrating an extension 
of the ground-state phase diagram, Fig. |31 to finite T, supported by theoretical 
calculations in the narrow- resonance regime |68| . The red dot is a tricritical point. 
For —-j^^ > —2.37, to the right of the QTP, the phase boundaries Pci and Pc2 extend 
to finite T, merging at a tricritical point. In the regime where the FFLO can be 
stable, one expects a sliver of FFLO phase along the Pc2 boundary. The schematic 
phase diagram depicts a putative regime of detuning, where at T 0, Pci > and 
Pc2 < 1. Such a regime does not appear within the mean- field approximation, but is 
found in QMC calculations [55] . 

its BEC-BCS crossover. While the simple BEC-BCS variational mean-field wavefunction 
gives a reasonable approximation for T = properties across the phase diagram, for 
the gas at finite temperature one must include the effect of thermally excited molecular 
pairs and fermionic single-particle excitations to obtain a reasonable T^. Because these 
excitations are expected to be strongly interacting around the unitary regime (FR point) 
a quantitatively trustworthy description remains a challenging open problem for both 
balanced and imbalanced gases. 

However, quantitative predictions for phase boundaries notwithstanding, it is likely 
that the form of the finite-T phase diagram, predicted by Ref. [68] via a controlled 
narrow resonance two-channel model analysis [271 EDI EHl [66] is qualitatively correct. 
The finite-T phase diagram, illustrated Fig. |5] is reminiscent of the phase diagram of 
^He-'^He mixtures. 

A number of new ingredients emerge at finite temperature. Because (in contrast to 
the fully-gapped BEC-BCS paired ground state) "spin"-ful quasi-particles can carry a 
finite imbalance, at finite temperature a chemical potential imbalance h always induces 
a nonzero species polarization. Thus at finite T the BEC-BCS singlet superfiuid state 
will no longer be distinguished from the SFm by a vanishing polarization. Therefore at 
nonzero temperature is a crossover, that at low temperature is expected to remain 
relatively sharp. On the BCS side of the resonance at low T the first-order transition out 
of the paired BCS superfiuid is expected to survive, though likely to reduce in strength. 
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Another finite temperature effect is tlie migration of tlie zero-temperature tricritical 
point of Fig. |3]into tlie P < 1 regime [681 170 1 1117] . as illustrated in Fig. |5l Finally, as we 
will discuss below, a strikingly more significant qualitative effect of thermal fluctuations 
is predicted to arise in the putative unidirectional FFLO states [T26] . The enhanced 
role of fluctuations in such periodically paired superfluids is associated with underlying 
rotational and translational invariance that is spontaneously broken. As we will discuss 
in more detail in the next section, fluctuations lead to universal power-law density 
correlations in the LO state, and a melting of the translational order (that can liberate 
fractional charge dislocations and vortices) and induce quantum nematic phases. 



3. Fulde-Ferrell-Larkin-Ovchinnikov superfiuid 

One of the more interesting alternatives to the resolution of the species-imbalance 
imposed frustration is the Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) [961 197] superfiuid. 
In this state the Fermi surface mismatch is accommodated by finite momentum pairing, 
that in its most generic form is a paired supersolid, that is somewhat of a holy grail 
in condensed matter physics. First proposed more than 45 years ago, this enigmatic 
paired state has been explored extensively in the context of solid state superconductors 
and high-energy physics systems [131 Hlj. However, with the exception of recent 
promising experiments on CeCoIup; [1181 1119] , the FFLO state has so far eluded definitive 
observations. 



3.1. Theoretical description 

The existence of the FFLO state is most easily detected, theoretically, through a 
transition from the polarized normal state upon reduction of the chemical potential 
difference, h (or equivalently, by reduction of the species imbalance). Assuming 
(and verifying a posteriori) a continuous (or weakly first-order) N-FFLO transition, 
a perturbative treatment via a Landau expansion in the pairing order parameter Aq, 
defined by A(x) = Aq^c*'^"'^, is justified. Standard analysis gives Landau free 
energy density 

Ugl = Yl ^Q^Q^Q + E ^Qi,Q2,Q3,Q4 Aq Aq. A*Q3 Aq, + . . . , (25) 
Q Qi,Q2,Q3 
where translational invariance constrains Eq to be a function of the magnitude of Q 
only, Q4 = Qi — Q2 + Q3, Aq is a Fourier transform of A(x), and 



3n 



Q 1, v'Q'-4h' h , VFQ + 2h 

1 — In — 1 In 

27/i 2 ^Bcs "^fQ vpQ — 2h 



(26) 



obtained either by a direct Landau expansion [97] or via expansion of the Bogoliubov-de 
Gennes (BdG) equation[66] to quadratic order in Aq. Simple analysis shows that eq 
has a minimum at |Q| ~ 1.81ABcs/^'yF ~ 0.58/^o, infinitely degenerate with respect to 
the detailed nature of the FFLO state, thereby leaving the quartic interactions to lift 
this degeneracy, i.e., to select the set of Q's defining a specific type of the FFLO state. 



Imbalanced Feshbach-resonant Fermi gases 



19 



As first demonstrated by LO, just below hc2 it is the ±Q periodic (cosine) state that is 
selected. 

However, the above treatment cannot be justifiably extended far below hc2, down 
to the transition to the uniform fully-paired superfiuid state near hd- There, the order 
parameter is no longer small and therefore must be treated nonperturbatively. Namely, 
a complete functional dependence of 7{[A(x)] may need to be taken into account, and 
an arbitrary form of A(x) (arbitrary set of Fourier harmonics Aq) must be considered. 

The findings of LO [97] motivate one to focus on a unidirectional FFLO order 
characterized by a colinear set of the Q„. Under this restriction, FFLO states fall 
into two, LO and FF universality classes. The LO (FF) states are characterized and 
distinguished by breaking (preserving) translational and preserving (breaking) time- 
reversal symmetries. The low-energy properties of such states can be well captured 
with a single Q pair (LO) and a single Q wave vector (FF) approximations, as we 
describe below. However, the analysis of energetics of these states far below hc2 requires 
the aforementioned fully nonlinear analysis. 

In principle this can be accomplished through a BdG treatment. However, its 
analytical solution (outside Id) is only tractable for a single Q FF state. This approach 
predicts a narrow sliver of FF state near hc2, that closes down at —l/kpa^ ^ 0.46, as 
the strongly coupled regime is entered [591 ESj . The most direct interpretation of this 
finding is refiected in our T = phase diagram in Fig J3(59l [66] , and is quite pessimistic 
for the observation of any of the FFLO states due to its narrow range of stability and 
the trap induced cloud inhomogeneity. 

However, an alternative, more optimistic scenario is suggested by early numerical 
treatments of 2d and 3d BdG equations |120l 1121] . exact Id solution |122j . density 
functional theory results |123] . and recent findings of a negative domain- wall energy 
in an otherwise uniform singlet BCS superfiuid |124l 1125] . These fully nonlinear (in 
A(x)) findings support the idea that a LO-type periodic superfiuid state of domain- 
walls (corresponding to a large set of colinear wavevectors in A(x)) may be significantly 
more stable down to the vicinity of hd. The LO state can be thought of as a periodically 
ordered micro-phase separation between the normal and BCS states, that thus naturally 
replaces the macro-phase separation ubiquitously found in the BEC-BCS detuning- 
polarization phase diagram. In Fig. [6l we show a proposed phase diagram, assuming 
that the SF-LO phase transition is indeed continuous and occupies a wider regime of 
the phase diagram than the mean- field prediction (see Ref. |125j ). 

Despite the extended history of FFLO states, only recently has the role of quantum 
and thermal fiuctuations been understood [126]. This study developed Goldstone mode 
models for the LO and FF universality classes and analyzed the corresponding beyond- 
mean- field phenomenology. Focusing on the more stable LO state [971 1120[ 11211 1125[ 1124] , 
it was demonstrated that in contrast to the conventional (uniform) superfiuid and 
FF states, the LO superfiuid exhibits two Goldstone modes, 9±q, corresponding to 
two primary order parameters A-tQ (with phases for high harmonics locked to these 
two). Equivalently, low-energy fiuctuations are characterized by the superfiuid phase. 
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Figure 6. Proposed schematic pliase diagram of AN/N vs. l/(A;Fas), showing LO 
liquid crystal phases replacing phase-separated (PS) regime. The three points are the 
same as in Fig. 3, defining the limits of stability of the surrounding phases. The lower 
panel shows symmetry-allowed 3d transition scenarios as a function of temperature 
to the normal- nematic (N-Nm), normal- isotropic (N-I), normal-smectic (N-Sm2Q), 
and "charge" -4 superfluid- nematic (SF4-Nm) phases. Currently, no microscopic 
calculations exist that assess which of the three possible scenarios are actually realized 
in imbalanced Fermi gases and in what range of detuning. 



9sc = {0+ + 6'-)/2 and the domain-wall phonon, u = Ogm/Q = (^+ — ^-)/2(5, with the 
LO order parameter (in a single ±Q approximation) given by 

Az.o(x) = 2AQe^^-(^) cos (Q ■ x + ^^,,„(x)). (27) 

The LO order parameter is a product of a superfluid order parameter and a 
unidirectional spontaneously oriented (along Q) Cooper-pair density wave, i.e., it 
simultaneously exhibits the off-diagonal long-range order and smectic order. 

Symmetry arguments supported by detailed microscopic calculations |126j show 
that the low-energy modes of the LO state involve two coupled smectics (with the layer 
normal taken to lie along xy), with moduli derivable from the BCS theory: 



n 



LO 



a form that is familiar from studies of conventional smectic liquid crystals, with 
rotational invariance encoded through the vanishing of the (V_lm)^ modulus and the 
specific form of the nonlinear elastic terms. The LO state is found to be a highly 
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anisotropic superfluid with the ratio of superfluid stiffnesses given by 

pj/pl = ^ (Aq/Abcs)' ^ Hh,2/h) « 1, (29) 

vanishing for h — )■ h~2. The FF state is found to be even more exotic, characterized 
by an identically vanishing transverse superfluid stiffness, a reflection of the rotational 
invariance of the spontaneous current to an energy-equivalent ground state. 

Thus, a resonant imbalanced Fermi gas conflned to an isotropic trap is a natural 
realization of a quantum (superfluid) liquid crystal that, unlike the solid state analogs, is 
not plagued by the underlying lattice potential that explicitly breaks continuous spatial 
symmetries. 

It has been demonstrated |126] that in 3d the long-range LO order is stable to 
quantum fluctuations, but is marginally unstable at any nonzero T. The resulting 
superfluid state is an algebraic phase, characterized by universal quasi-Bragg peaks 
and correlations that admit an asymptotically exact description. In contrast, putative 
crystalline LO phases [T27] with multiple noncolinear ordering wavevectors are stable 
against thermal fluctuations. The resulting state is also found [T26] to exhibit an unusual 
topological excitation - a half vortex bound to a half dislocation - allowed by the above 
form of A (x) in Eq. ([27]). 

In 2d, at nonzero T, the state is also an algebraic phase, exhibiting universal 
power-law phonon correlations, controlled by a nontrivial exactly calculable fixed point. 
However, it displays short-range positional order with Lorentzian structure-function 
peaks, and is thus unstable to the proliferation of dislocations. The resulting state is 
either a "charge" -4 (four-fermion condensate) superfluid (with order parameter A^Jc ~ 
A^ ^ iAgC^^^^'^), or a non-superfluid nematic, depending on the relative energetics 
of the aforementioned integer and half-integer vortex-dislocation defects. The latter 
normal nematic state is a (complementarity described |128] ) deformed Fermi surface 
state |115|, 1116] , that is another exotic candidate for a strongly interacting, imbalanced 
Fermi gas. 

The unbinding of various combinations of topological defects (dislocations and 
vortices) yields predictions for a number of other interesting phases [T26] . These include 
two types of aforementioned orientationally-ordered but statistically homogeneous 
(dislocations are unbound) nematic states, one where superfluidity is completely 
destroyed ("normal" nematic, N-Nm) and another in which 4-atom superfluid order 
accompanies the orientational nematic order ("superfluid" nematic, SF4-Nm). Another 
state is the "normal" smectic phase, N-Sm2Q, emerging from the LO state by 
proliferation of vortices, that destroy the superconducting phase coherence, but retain 
the periodic density layering. The three possible symmetry-allowed scenarios for 
transitions between these liquid crystal phases as a function of temperature are 
illustrated in the lower panel of Fig. [61 Currently, no microscopic calculations exist 
that assess which of the three possible scenarios are actually realized in imbalanced 
Fermi gases and in what range of detuning. 

A complete low-energy description of these LO liquid crystal states must take the 
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gapless fermionic excitations into account, in addition to the above bosonic Goldstone 
modes. As with other analogous problems [115 j, these are expected to lead to the 
Landau-like damping of the Goldstone modes 6sc, u, and a finite fermionic quasi-particle 
lifetime. So far, this difficult problem has not been addressed. 



3.2. Experimental predictions 

The possibility of the observation of conventional FFLO states has been discussed at 
length [621 1129i [59| [66] . Bragg peaks (refiecting spontaneous periodicity) observed in 
time of fiight [66] experiments and the spectra of collective breathing modes [T30] are 
two prominent signatures. It has been recently emphasized that the FFLO states may 
be energetically stabilized in a quasi-ld geometry |131] . 

In the inhomogeneous environment of the trap there are a number of constraints on 
the observation of FFLO states and the fiuctuation phenomena discussed above. Near 
hc2 the LO period Xq = 271 /Qq is bounded by the coherence length (that near unitarity 
can be as short as ~ R/N^^^, where R is the trapped condensate radius and is the 
total number of atoms), and thus -C R and in this regime the trap can be treated via 
a local density approximation (LDA). For Xq -C R, LDA predicts a weak pinning of 
the LO smectic, that can be estimated via finite size scaling, with trap size R cutting 
off (m^) ~ rjlog^R/ Xq), leading to (A^o) ~ {^q/R)^ ^ 1 that no longer truly vanishes 
although it is still strongly suppressed. We expect the predicted strong-fiuctuation 
effects to be experimentally accessible. We note, for example, that Kosterlitz-Thouless 
phase fiuctuation physics has been reported in 2d trapped superfiuids[lT], despite the 
finite trap size. However, a more detailed analysis of the trap effects, necessary near hd 
for a quantitative comparison with experiments, remains to be done. 

Other experimental signatures of the LO state include vortex fractionalization, 
where the basic superfluid vortex is half the strength of a vortex in a regular paired 
condensate. This should be observable via a doubling of the vortex density in a rotated 
state. Also under rotation, the high superfluid anisotropy Eq. ( I29l) is expected to lead 
to an imbalance-tunable strongly anisotropic vortex core and a lattice that is highly 
stretched along Q. Bragg peaks in the time-of-flight images can distinguish the periodic 
SF2-SmQ (superfluid smectic) state from the homogeneous SF4-Nm (superfluid nematic), 
which are in turn distinguished from the N-Sm2Q and N-Nm (normal smectic and 
nematic) by their superfluid properties, periodicity, collective modes, quantized vortices, 
and condensate peaks. Thermodynamic signatures can identify the corresponding phase 
transitions. 



4. Local density approximation 

One of the principal experimental probes of ultra-cold atomic gases is the measurement 
of the local atom density of the trapped cloud. Thus, quite generally, detailed theoretical 
predictions for the phases and phenomena of cold polarized Fermi gases must account 
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Figure 7. Schematic depiction of the shell structure of polarized spherical cold-atom 
clouds, with the phases occurring with increasing radius in a trapped polarized Fermi 
gas, within the LDA, (a) in the BCS and unitary regimes and in (b) the BEC regime 
where the SFm phase exists. 



for the effect of the trapping potential Vrir). In recent years, there have been numerous 
theoretical investigations of the density profiles of trapped polarized fermion gases (see, 
e.g., Refs. |59l [1321 [1331 IIHI [1351 [1361 here we 

shall focus on the predictions of the simplest method for handling the trapping potential, 
namely the local density approximation (LDA). 

Much like the WKB approximation, the LDA corresponds to using expressions for 
the bulk system, but with an effective local chemical potential fi{r) = fi — Vrir) in 
place of n- Systems for which the LDA holds, then, are of particular interest as they 
are relevant for the comparison to correlated condensed matter. The validity of the 
LDA relies on the smoothness of the trap potential, namely that Vr(r) varies slowly 
on the scale of the longest physical length A (the Fermi wavelength, scattering length, 
effective range, etc.) in the problem, i.e., {X/VT{r))dVT{r) /dr <^ 1. Its accuracy can be 
equivalently controlled by a small parameter that is the ratio of the single particle trap 
level spacing 6E to the smallest characteristic energy of the studied phenomenon 
(e.g, the chemical potential, condensation energy, etc.), by requiring 6E/Ec ^ 1. 

Because of its more direct current experimental relevance, in this section we focus 
on the single-channel model, Eq. ( [TT]) . The generalization of this model to a trap is 
straightforward: 

H=J2[ d'rciiv) [ - ^ + (VT(r) - /i.)] c,(r) + A /" rfVcjcjc/^, (30) 

where Co-(r) is a fermionic field operator with Fourier transform Cko-- Henceforth, to 
be concrete, we shall focus on an isotropic harmonic trap Vr(r) = Vrir) = |mfi^r^, 
although this simplification can be easily relaxed to handle an arbitrary anisotropic 
trap. The LDA assumes that, locally, the system is taken to be well approximated as 
uniform, but with a local chemical potential given by 

/x(r) = fi- ^mfllr'^, (31) 
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Figure 8. Plot of the local fermion densities, n-f- (dashed) and (solid) , as a function 
of radius and normalized to the overall density uq, as a function of position normalized 
to the P = Thomas-Fermi radius Rq, in the BCS regime with (feplasl)^^ — 1-5, from 
Ref. [55]. 



where the constant /i is the true chemical potential (a Lagrange multiplier) still enforcing 
the total atom number N. The spatially-varying spin-up and spin-down local chemical 
potentials are then given by 

/it(r) = /i(r) + h, (32) 

f^iir) = Kr) - h, (33) 

with the chemical potential difference h uniform. Thus, within the LDA we approximate 
the energy density by that of a uniform system with spatial dependence (via the trap) 
entering only through /i(r). The ground state energy is then simply a volume integral of 
this energy density. Within LDA, the phase behavior as a function of chemical potential, 
/i, translates into a spatial cloud profile through /i(r), with critical phase boundaries fic 
corresponding to critical radii defined by fic = /^('^c, h). As we first predicted [59], this 
leads to a shell-like cloud structure (illustrated in Fig. [7]) that has subsequently been 
observed experimentally [551 ESI EZ] . 

Below we briefly review these shell structures using the LDA, following our earlier 
work [66] • We note, however, that throughout our discussion, sharp (discontinuous) 
features that arise are an artifact of LDA and are expected to be smoothed on 
microscopic (Fermi wave-) length scales by the kinetic energy (or, surface tension |137j ). 

We begin with the analysis in the BCS regime, where /i > 0. Qualitatively, the 
unitary regime will behave analogously to the BCS regime. Neglecting the FFLO phase, 
the fixed /x and h phase diagram has a first-order phase transition to the normal phase 
at 

he ^ --^ABEC-Bcs(/i), (34) 

^ 4y2e- V(r) exp [- ^^^y^l , (35) 

L47rasA^(/i(r))J 

separating a polarized normal phase (when h > he) and a paired superfiuid phase 
(when h < he). Since ABEC-Bcs(Ai) is an increasing function of yii, this leads to the LDA 
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Figure 9. Phase diagram, in the grand-canonical ensemble at fixed /i and h, with h 
normalized to the molecular binding energy Ei, = and — 2fi + Ei, the effective 
molecular chemical potential. The thick grey line is a first-order transition, and the 
thin black lines are continuous transitions. The red point is a tricritical point and the 
purple point is a critical end point. 



prediction of a fully paired superfluid core at small r (where /i is larger) surrounded by a 
population imbalanced normal phase, as illustrated in Fig. [7^. This yields a remarkably 
simple behavior for the spin-'l' and spin- J, densities n^^i, plotted in Fig. |8l with = 
in the central core (the fully-paired superfluid phase), and ^ rii in the outer shell 
(the imbalanced unpaired Fermi liquid phase). We note here that including the FFLO 
phase, that occurs near the SF-N phase transition, would imply a thin shell of FFLO 
within the LDA yielding oscillations of the local pair amplitude near the SF-N interface, 
as also found in recent beyond-LDA analysis based on the BdG equations |133j . 

Turning to the BEG regime. Fig. |9] shows the homogeneous-case phase diagram at 
fixed n and h that is required to construct the phases occurring with increasing radius 
within the LDA. The structure of this phase diagram is somewhat nonintuitive; this is 
because at fixed densities the chemical potential in the EEC regime changes rapidly with 
external parameters (like the detuning), as also seen in the balanced case, see Fig. [T|d. 
Nonetheless, the resulting shell structures are the natural generalization of the BCS 
regime results: Within the LDA, possible sequences of phases with increasing radius 
correspond to constant h trajectories, with decreasing /i, in the phase diagram of Fig. |9l 
Thus, we see that the possible sequences of phases with increasing r are: pure SF, SF— >■ 
N, or SF— j-SFm— N, with the latter case illustrated in Fig. [3d. 

The characteristic LDA profiles in the BEG regime for the local pair density nm{f) 
and magnetization M(r) are shown in Fig. [10] for the sequence of phases SF— j-SFm— >■ N. 
Thus, as illustrated in Fig. [TJ in contrast to the BCS regime here the cloud can exhibit 
a thin shell of the SFm phase with finite polarization and a pairing gap. The three 
characteristic radii of importance are: (the SF— j-SFm transition), Rtf (the Thomas- 
Fermi radius, where the SFm—)- N transition occurs) and i?/2, the outer boundary 
of the cloud. As we discuss in the next section, the LDA has been used by the 
Ketterle group |108] to infer the phase diagram for a homogeneous system. However, 
interesting physics also holds for systems in which LDA is violated. As discussed 
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Figure 10. Typical pair density, rim oc |Ap, and local magnetization, M, each 
normalized to a characteristic density scale EbN^Ei,), with N{E]y) the density of states 
measured at the binding energy, from Ref. [SS] in the BEC regime. 



by DeSilva and Mueller |135] . the early Rice experiments [56] showed axial density 
profiles (corresponding to the density, integrated over two directions, as a function of 
the coordinate z along the long axis of the cloud) that apparently violated the LDA. 
Indeed, the violation of the LDA can be directly seen in the data shown in Fig. 1 of 
Ref. [56], since the fermion densities do not follow the contours of the large aspect 
ratio trap. Additionally DeSilva and Mueller showed (see also Ref. [136 J) that for 
LDA-type cloud shapes of the form of Fig. [TJ the axial magnetization must decrease 
monotonically with increasing z. A principal culprit in the breakdown of LDA is the 
surface tension associated with the SF-N interface |137] . This issue was studied in detail 
in Ref. |143] in an attempt to understand the density profiles from the Rice group [56] 158]. 
However, it was found theoretically that although surface tension indeed deforms the 
minority cloud in a way qualitatively consistent with the experiments, the size of the 
deformation is too small to account for experimental observations based on microscopic 
estimates of surface tension magnitude. This leaves a quantitative understanding of 
the Rice experiments an open question; one possible important ingredient that needs to 
be included are fluctuations around mean-field theory (known to be important in the 
unitary regime), as discussed in recent work |142j . 

5. Recent experiments 

We now turn to a more comprehensive discussion of recent experimental developments 
in imbalanced Feshbach-resonant fermionic atomic gases. These have primarily been led 
by the Ketterle group at MIT and the Hulet group at Rice, launched by two studies of 
^Li that appeared in Science in 2006 [55l [56] . 

The MIT experiments [55] reported a destruction of superfluidity (probed by the 
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appearence of vortices in response to an imposed rotation) and condensation (probed 
by a fit to the cloud shape after sweeping to the BEC side of the resonance) with 
increasing imbalance. As shown in Fig. 1 of Ref. [55], a rotating polarized Fermi gas 
exhibits a distinctive inner superfluid core (defined by a region to which the induced 
vortex lattice is confined) surrounded by a vortex-free shell of imbalanced normal cloud. 
Consistent with this identification, Fig. 2 of Ref. [55] shows that the number of vortices 
A^^ decreases with increasing imbalance for different values of the detuning. On the 
basis of the standard relation |144] Ny = 2mQA/TTh, we thus observe that simply 
measures the cross-sectional area A (normal to the axis of rotation) of the superfluid 
core as a function of the rotation rate Q. Although qualitatively this data (the shrinking 
of the superfluid core with increasing imbalance) is consistent with current theoretical 
predictions, these observations remain to be quantitatively understood. Zwierlein, et 
al. [55] also reported evidence of phase separation, as reflected in a suppression of the 
local polarization (species imbalance) near the center of the cloud, a signature of a core 
that is fully-paired into a BEC-BCS singlet superfluid. These were confirmed by in situ 
measurements of individual species density profiles using an exquisite phase-contrast 
imaging and 3d image reconstruction technique |57]. Based on these observations they 
constructed a phase diagram [55] with a phase boundary separating the normal phase 
and superfluid phases (the latter identified with normal-superfluid phase separation), in 
qualitative agreement with the mean-field theoretical phase diagram [591 166] . 

The initial Rice experiments [56] studied the imbalanced ^Li Fermi gas by comparing 
the density profiles of the two fermion species. Consistent with the MIT experiments 
they found that, with increasing species number imbalance, the imposed imbalance is 
pushed to the outside shell of the cloud. Furthermore, by measuring the dependence 
of the radii of the Thomas-Fermi clouds for the two fermion species and showing that 
the ratio of the majority to minority radius, R^/R^, increases with increasing imbalance 
P, Partridge, et al. showed evidence for phase separation in the unitary regime. From 
these measurements, the Rice group found a putative nonzero lower-critical polarization 
for phase separation with Pd ~ 0.09, observing that for P < Pd the radii of the 
two species coincided, R^ = Ri, despite a nonzero imposed imbalance (see Fig. 3 of 
Ref. [56]). We are not aware of a ground state (within a comprehensive theoretical 
picture) that can accomodate this observation of a nonzero Pd, i.e., a state that is both 
fully paired and imbalanced; naive candidates like the homogeneous polarized superfluid 
SFm (appearing on the BEC side) [591 [66] and a putative p-wave paired state [145] would 
fail toward the edge of the cloud, where on general grounds a polarized normal state 
with distinct radii is expected. Furthermore, a reasonable understanding of this data 
was obtained by Chevy [81] within a simple unitary model of phase separation (i.e., 
with distinct majority-minority radii) between a fully-paired unpolarized superfluid and 
a fully-polarized normal state, using universal forms for the equations of state for these 
phases. Taking these together with the fact that no such finite Pd was ever observed by 
the MIT group [57] suggests that the putative (small) nonzero Pd in Rice measurements 
may be an experimental artifact. An alternative possibility is that it is a consequence 
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of a nonzero temperature T ^ O.lTp, that on one hand allows a finite polarization 
(through thermally activated quasi-particles) and on the other entropically suppresses 
phase separation, thereby allowing a nonzero Pci- One support for the latter scenario 
is the fact that in the later Rice experiments ^8] a vanishing P^i was found at reduced 
temperature, T ^ O.OSTp. Detailed experiments further clarifying this issue would be 
highly desirable. 

A number of subsequent experiments further explored the onset of phase separation 
in such imbalanced Fermi gases. Zwierlein, et al. |55] studied the transition between 
the unpaired and superfiuid phases by observing column density profiles in imbalanced 
Fermi gases, finding the upper-critical polarization of Pc2,trap — 0.70 at the unitary point 
for the trapped cloud. Due to the presence of the trap, this value cannot be directly 
compared to the upper-critical polarization Pc2 for the bulk, uniform system, the latter 
estimated within mean-field analysis to be P^^^ ^ 0.93. However, within the LDA, quite 
generally 



This bound that can be understood by considering a trapped imbalanced Fermi gas 
at Pc2,trap = M/n = AN/N, with magnetization M = AN/V {V the system's 
volume) and atom density n = N/V. Within the LDA, at Pc2,trap, by definition the 
fully-paired inner superfiuid core has shrunk to zero radius. Thus, at the center of 
the trap the local polarization is given by the bulk upper-critical polarization, i.e., 
P{r = 0) = M(0)/?T,(0) = Pc2; where M(r) is the local magnetization and n{r) is the 
local density. Now, since /i(r) = /i — \mVL^r'^ decreases with increasing radius, and the 
chemical potential difference h is r-independent, M(r) and n(r) respectively increase 
and decrease with increasing radius, implying 



and giving the bound in Eq. ( l36|) . Thus, the observation of Pc2,trap — 0.70 in Ref. |55] is 
a measure of a lower-bound on the inaccuracy of the mean- field estimate for Pc2- 

The imaging technique in such experiments naturally integrates the local atomic 
density along the axis perpendicular to the camera's imaging plane. To deconvolve this 
projection. Shin, et al. [57] used the inverse Abel transform [22] to extract the local 
three-dimensional atom densities. This development allowed a clear identification of the 
paired superfiuid core (within which the local imbalance identically vanishes, at T = 0) 
in the phase separation regime, and the measure of the upper-critical polarization as 
a function of detuning, finding Pc2,trap = 0.77 at unitarity, close to unity in the BEG 
regime, and decreasing in the BCS regime. 

One of the most significant differences between the MIT and Rice observations is 
the shape of the interface between the balanced superfiuid core and the imbalanced 
outside shell. While the MIT group found an interface whose shape mimics the trap 



Pc2,trap > Pc2- 



(36) 




(37) 



(38) 
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(see, e.g., Fig. 6 of Ref. [SZ]), the Rice group observed that the imposed imbalance 
is accomodated by an enhanced local polarization in the axial poles of the cloud |56j . 
thereby at low T deviating from the trap's highly anisotropic shape. This discrepancy, 
which to date remains to be quantitatively understood, is believed to be related to 
the very different trap aspect ratios in the MIT {fr/fz = 5.6 in Ref. [5^) and Rice 
[fr/ fz = 49 in Ref. [56]) experiments, the latter exhibiting a nearly one-dimensional trap 
geometry. This issue was studied in more detail in Ref. [58], where the distorted shape 
(attributable to a nonzero surface tension |135] associated with the first-order transition 
and corresponding phase separation) of the paired inner core was confirmed, and the 
dependence of the aspect ratios of the minority and majority species fermionic clouds as a 
function of the imposed imbalance was measured. These later Rice experiments [58] also 
found a surprisingly large value of the in-trap upper critical polarization, Pc2,trap ~ 0.9. 
Subsequent Rice experiments jl46] explored the effects of trap aspect ratio, and found 
that the LDA-violating deformations decrease with decreasing aspect ratio, with no 
deformations observed for aspect ratios below 6. Furthermore, the critical polarization 
was found to be approximately P ~ 0.75, consistent with the MIT measurements and 
Monte Carlo calculations. 

However, recent Paris experiments [71] studied the phase diagram of an imbalanced 
Fermi gas with around 10^ atoms and trap aspect ratio of around 20, conditions in 
between the Rice and MIT systems. In contrast to the above appealing resolution of 
the discrepancy between the Rice and MIT observations (due to a very different aspect 
ratio), these Paris experiments saw no appearance of LDA-violating features or surface 
tension effects. Furthermore their results are in a very good agreement with MIT's, both 
qualitatively (observing a fiat-top distribution of the density difference, that can arise 
from a fully-paired superfluid core within the LDA), and quantitatively (with respect to 
the critical polarization Pc2)- Thus, a resolution of these puzzling discrepancies remains 
an open issue. 

The observed matching of the majority and minority clouds' aspect ratios to the 
trap's aspect ratio in the MIT experiments indicates that, for these experiments, the 
LDA is well satisfied. The MIT group then took advantage of the LDA to extract 
the T = and finite T phase diagrams \108\ 1109] for a bulk system, facilitating a 
more direct comparison with theoretical studies of a homogeneous gas. Focusing on 
the unitary limit, these experiments found a finite-T phase diagram of the standard 
tricritical form (i.e., like Fig. [5l resembling other tricritical systems, such as ^He-^He 
mixtures |104j . metamagnets |147j . and thin- film superconductors |148j ). 

As illustrated in Fig. |5l a tricritical point separates a second-order (at high T) 
from a first-order (at low T) transition that, in the present case, is between a paired 
superfluid and a polarized normal state, with the latter extending to a regime of phase 
separation at fixed polarization. This general structure is well captured by mean-field 
theory for the narrow- resonance model [68], however the shape of the phase diagram 
in the experimentally-relevant wide-resonance regime remains an open question. Finite 
temperature analysis beyond mean-field theory |149] finds a more complicated phase 
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diagram topology, with the regime of phase separation persisting for temperatures higher 
than the tricritical point. 

The principal experimental evidence for the expected tricritical point comes 
from the observations [58l 1108] that at higher temperatures the breakdown of LDA 
(as signaled by the distinct cloud and trap shape anisotropics) ceases. This was 
reasonably interpreted as the disappearance of the surface tension (associated with 
phase coexistence at first-order transition) as the low-T first-order superfluid-normal 
transition converts to a continuous one for T above the tricritical temperature. This 
is further supported by observations that at these higher temperatures (above the 
tricritical temperature) the expected Clogston hmit of Pc2,trap ~ 0.75 is recovered. A 
more detailed check on the consistency of the interpretation can come from a comparison 
with the predictions for the cloud shape, heat capacity, and local polarization controlled 
by the mean-field tricritical point |117] . which are predicted to follow distinct power 
laws [T501IT5T] . 

As mentioned above, in later work |109] . Shin, et al. also used the LDA to extract 
the T = phase diagram of an imbalanced Fermi gas. The resulting experimental phase 
diagram agrees qualitatively with the theoretical one displayed in Fig. [3] [591 ES] , though 
no experimental evidence for the putative FFLO phase exists to date, nor has the low-T 
SFjvi-to-PS phase boundary (Pci) been mapped out. The observed value of Pc2 — 0.36 
at unitarity (reported in Ref. |108] . with the phase diagram of Ref. |109] showing a 
slightly lower value) is consistent with the above bound Pc2 < -Pc2,trap and points to 
a large quantitative error in the mean-field estimate (although including the leading 
correction [83] yields a more accurate estimate, Pc2 — 0.302) This measurement 
is consistent with the quantum Monte Carlo based prediction by Lobo, et al |112] . who 
find n-^/n^ = 0.44 at the transition, corresponding to Pc2 = 0.38 and within LDA giving 
a universal value of Pc2,trap = 0.77, in agreement with the MIT [57] and latest Rice |146] 
measurements. 

One remarkable aspect of the experiments in Ref. |109] is that the density profiles of 
the trapped imbalanced Fermi gases were quantitatively captured within a simple boson- 
fermion mixture model (with bosons representing the tightly bound molecular pairs 
and fermions representing the excess spin- up fermions), using the vacuum molecule- 
fermion |152] (at/ = 1.18a) and molecule-molecule |153[ 185] {ubb = 0.6a) scattering 
lengths. Since the latter values are only valid in the dilute deep-BEC regime, it 
is surprising that such an accurate fit to the data held close to unitarity. It will 
be interesting, in future work, to test the limits of the boson-fermion model used in 
Ref. |109] and to look for deviations from this simple picture. 

Before turning to the last topic of this section, which is the development of RF 
spectroscopy probes of imbalanced Feshbach-resonant Fermi gases and the strongly- 
imbalanced polaron regime, we briefly comment on recent Rice experiments |154] . that 
have attained the quasi one- dimensional conditions for which the Id FFLO state is 
expected. These experiments utilize a two-dimensional optical lattice to create an 
array of weakly-coupled, large aspect ratio tubes (with Ur/oJz = 1000), that are 



Imbalanced Feshbach-resonant Fermi gases 



31 



therefore effectively in a one- dimensional limit. Previous theoretical work [131] and 
general arguments suggest that in this limit the FFLO state should robustly appear. 
Moreover, uniform one-dimensional imbalanced Fermi gases can be solved exactly via 
Bethe ansatz |155[ I156j . supporting these claims. Liao et al |154] found that, combining 
such Bethe ansatz results with LDA (to account for the spatial dependence along the 
tube quantitative agreement with their trapped density profiles was obtained. 

Consistent with earlier theoretical predictions |131] . these experiments found that in 
contrast to the 3d case, the superfluidity is weakest and the imbalance (magnetization) 
is largest at the center of the trap where the chemical potential is smaller. This is a 
consequence of the reduction in the Id density of states with the increasing chemical 
potential (contrasting with the opposite behavior in 3d), resulting in the appearance of 
the balanced superfluid at small fi, and the imbalanced superfluid (the Id equivalent of 
the FFLO state) at large fi at the center of the trap. An important goal for future work 
is a direct observation of the expected modulation of the FFLO state. This should be 
reflected in quantities such as the local pairing correlations, as found in recent QMC [157] 
and density-matrix renormalization group |158[ 1159] studies of one-dimensional systems. 
Finite momentum Bragg BEG peaks (akin to those appearing in a condensates in an 
optical lattice, but here without an imposed periodic potential) are also expected to 
characterize the time-of-fiight measurements of the FFLO states. 

We now discuss RF spectroscopy [H], which has been fruitfully used to probe 
resonant imbalanced Fermi gases, with the particular focus on the strongly interacting 
normal state driven by a high imbalance and/or increased temperature [72l 1160] . RF 
spectroscopy studies of strongly imbalanced Fermi gases have a close condensed-matter 
counterpart, namely tunneling in the paramagnetic phase above the Clogston limit in 
thin-film superconductors |16H 1162] . Instead of removing or adding electrons through 
an insulating barrier, in RF spectroscopy the RF field induces transitions from one 
of the two interacting hyperfine species of the atomic cloud to a third unpopulated 
hyperfine state, and this rate as a function of the RF frequency u is measured by 
detecting the number N^lu) of the third species. In the simplest interpretation the 
signal N^i^uj) measures the spectral function of the species undergoing a transition to 
the third hyperfine state. 

However, as was recently emphasized [90l I163[ I164[ I165[ I166[ I167[ 1168] . this 
interpretation is complicated by final state interactions (a crucial issue that has received 
much recent theoretical attention [901 11631 [1641 11651 11661 1167[ 1168] ) of the third hyperfine 
state that are strong in ^Li [72| 1160] , and a significantly more involved theoretical 
analysis is required. (In contrast, in final state interactions are weak and a direct 
comparison with the fermionic spectral function is possible |169] ). 

The key features of the RF signal are the width and the shift of the peak in iV3(ci;) 
relative to that in a cloud of noninteracting atoms. Based on earlier studies [T7] the 
observed peak shift was originally associated with the pairing gap. An observation of a 
such shift for a strongly imbalanced gas in the regime (T,P,/ci?a) of the phase diagram, 
where the gas did not display a condensate peak, was taken as evidence for an exotic 
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paired non-superfluid state [72j. This is reminiscent of photoemission in the cuprate 
high-temperature superconductors, which exhibit a pseudogap |15], loosely taken as 
evidence of local pairing correlations |170[ I17H I172[ [38] in the non-superconducting 
state above T^. However, a detailed analysis [90] (including final state interactions 
and controlled by a systematic l/./V/-expansion[82l [83] ) of the imbalanced resonant gas 
demonstrated that the observed RF shift and widths can be qualitatively understood 
in terms of a conventional Fermi-liquid state, albeit a strongly renormalized one. The 
latter is characterized by large self-energy corrections leading to a modified Migdal 
discontinuity in the momentum distribution function (single-particle residue), Z < 1, 
the enhanced effective mass, m* > m, and a large chemical potential shift responsible for 
the observed shift in the RF peakjOO]. Subsequent experiments from the same group [T60] 
have successfully reinterpreted the original observations in terms of this theoretically 
advocated Fermi-liquid picture, though a detailed quantitative description of the data 
is still lacking. 

Finally, as mentioned earlier recent experimental studies [731 have provided 
additional insight on this interesting system in the strongly imbalanced limit of P — ?■ 1, 
corresponding to a single minority (e.g., spin down) atom moving and resonantly 
interacting with a fully polarized Fermi sea of majority atoms (spin up). Schirotzek, et 
al [73] find evidence for the predicted transition [173[ 1113] from a polaron-like state of 
the minority atom interacting with the majority Fermi sea, characterized by a nonzero 
fermionic pole with a residue Z < 1, to a diatomic bosonic molecule in the Fermi 
sea, with Z — )■ 0. Subsequent work by Nascimbene, et al. [71] studied the collective 
oscillations of the strongly imbalanced Fermi gas, extracting the effective polaron mass 
m* ~ 1.17m. These measurements are in good agreement with simple variational 
calculation based on the Chevy ansatz [SH 1111] and diagrammatic Monte- Carlo [T73], 
that both give m* = 1.17m, and another recent calculation [111] that yields m* = 1.20m. 

6. Concluding remarks and future outlook 

In this brief review, we have discussed recent developments in the study of imbalanced 
resonant Fermi gases. As we described, their phenomenology is summarized by the phase 
diagram as a function of imbalance (polarization P) and interaction [l/kpa, controlled 
by FR detuning). The zero-polarization paired-superfiuid ground state, exhibiting the 
well-studied BEC-BCS crossover, is readily destablized by an imposition of a finite 
species imbalance, on the BCS side, to a regime of phase separation between the paired 
superfluid and a polarized normal state. In contrast, sufficiently deep in the EEC regime 
side this is preempted by a polarized homogeneous superfluid, consisting of molecular 
Bose-condensed pairs and a fully polarized Fermi sea of the remaining majority- species 
atoms. Beyond the so-called Clogston limit of the upper-critical imbalance, the strongly 
interacting normal state is obtained. The properties of these phases and regimes have 
been vigorously studied through measurements of the number and imbalance density 
profiles, collective modes, as well as the RF spectroscopy as a function of interaction 
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strength (controlled by FR detuning), imbalance and temperature, finding qualitative 
and at times quantitative agreement with theoretical predictions. 

Despite this coherent picture that has emerged, there are number of interesting 
issues that remain to be explored. Firstly, more detailed and quantitative studies 
of the phase boundaries are needed, in particular a determination of the SF^ phase 
boundary P^i and the upper and lower critical polarizations as a function of detuning 
and temperature. In particular, it is not known whether the finite-T phase diagram has 
the standard tricritical form for all detunings. 

Secondly, more detailed investigations of phases are needed. As discussed in the 
manuscript, the polarized superfluid SFm, appearing in the BEC regime at T = 0, is 
described, at the simplest level, by a simple Bose-Fermi mixture of molecular pairs and 
single-species Fermi sea. It would be interesting to explore further the extent to which 
this picture is only an approximation and what the limit of its validity is, particularly 
as the system is taken closer to the unitary point, where s-wave scattering length grows 
and the molecules can no longer be treated as point bosons. 

There is also much that remains to be elucidated about the nature of the strongly 
interacting normal state at high polarization. Undoubtedly, RF spectroscopy and its 
recently-developed momentum-resolved extension |169] will be central to detailing its 
properties as a function of detuning and temperature. 

But probably the biggest remaining enigma is the long sought-after FFLO state, or 
more generally an imbalance-driven paired superfluid state that spontaneously breaks 
translational and orientational symmetries. On the theoretical side, the form of the 
periodic state (unidirectional or crystalhne and of what type |127] ). the range of stability 
of the state in the detuning-polarization-temperature phase diagram, as well as the 
nature of the transition to it from the paired BEC-BCS superfluid remain as important 
open questions. Furthermore, as recently demonstrated the simplest LO state admits 
fractional topological defects (a 1/2- vortex bound to a 1/2-dislocation), and exhibits 
enhanced fluctuations that lead to a variety of other quantum liquid crystal phases and 
other rich phenomena [126]. 

As discussed above, reliable analytical treatments around /ic2|SZ] and those 
specializing to the FF (single Q) state[96], and more recent ones that extend these 
beyond the deep BCS regime [591 166] suggest that the FFLO-type states are confined 
to a narrow sliver near hc2, just below the normal state. Whether this is generic or an 
artifact of current approximations remains an open question. 

One of the reliable limits that sheds some light on this question is the case of one 
dimension, where there is an exact solution of the Bogoliubov-deGennes equations, 
predicting a LO type ground state over a large portion of the phase diagram [120[ 
I155[ 1156] . This one- dimensional regime is also particularly amenable to numerical 
treatments [1571 11581 11591 11741 1175] and exact Bethe-ansatz methods [1551 1156] . In fact 
in one dimension, on general grounds the LO type state is generic for a finite imbalance, 
indistinguishable from the Id SFm ground state. Thus, the existence of a FFLO type 
state in one dimension is theoretically indisputable. Stimulated by this and supported 
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by quasi- Id mean- field analysis |131] . recent efforts have focused on a realization of the 
FFLO state in an array of weakly coupled Id imbalanced resonant Fermi gases, created 
by a strong 2d optical standing wave potential |146] . In the time of flight, this state is 
expected to exhibit spontaneous Bragg peaks associated with the simultaneous presence 
of superfiuid coherence and periodic positional order. Undoubtedly, many unanticipated 
surprises will emerge as a result of these and many other future studies. 
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